Phase Transitions in Small Systems: 
Micro canonical vs. Canonical Ensembles 



Jorn Dunkel 

Institute for Physics, Universitdt Augsburg, Universitdtsstrafie 1, D-86135 

Augsburg, Germany 

Stefan Hilbert 

Institute for Physics, Humboldt- Universitdt zu Berlin, Newton- Strafie 15, D-124-89 

Berlin, Germany 



Abstract 

We compare phase transition(-like) phenomena in small model systems for both 
microcanonical and canonical ensembles. The model systems correspond to a few 
classical (non-quantum) point particles confined in a one-dimensional box and inter- 
acting via Lennard- Jones- type pair potentials. By means of these simple examples it 
can be shown already that the microcanonical thermodynamic functions of a small 
system may exhibit rich oscillatory behavior and, in particular, singularities (non- 
analyticities) separating different microscopic phases. These microscopic phases may 
be identified as different microphysical dissociation states of the small system. The 
microscopic oscillations of microcanonical thermodynamic quantities (e.g. temper- 
ature, heat capacity, or pressure) should in principle be observable in suitably de- 
signed evaporation/dissociation experiments (which must realize the physical pre- 
conditions of the microcanonical ensemble). By contrast, singular phase transitions 
cannot occur, if a small system is embedded into an infinite heat bath (thermostat), 
corresponding to the canonical ensemble. For the simple model systems under con- 
sideration, it is nevertheless possible to identify a smooth canonical phase transition 
by studying the distribution of complex zeros of the canonical partition function. 
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1 Introduction 



One of the most intriguing thermodynamic properties of various macroscopic 
systems is their abihty to undergo phase transitions (PTs) if one or more 
control parameters pass certain critical values [1,2]. The first systematic clas- 
sification scheme for macroscopic PTs was proposed by P. and T. Ehrenfest [3] 
in 1912 already. After further pioneering work by Mayer et al. [4-7], Yang and 
Lee [8, 9] elucidated the mathematical essence underlying PTs in the grand- 
canonical ensemble by studying the distribution of complex zeros (DOZ) of 
the grandcanonical partition function. Later on, Fisher [10] and Grossmann et 
al. [11-13] employed a very similar approach to analyze the temperature de- 
pendence of phase transitions in the canonical ensemble. Recently, further sig- 
nificant progress in the understanding of critical phenomena has been achieved 
by studying the connection between PTs and phase (or configuration) space 
topology [14-20]. 

Formally, the seminal contributions [3-13] have in common that, in the spirit 
of traditional thermodynamics, they refer to macroscopically large systems; 
more exactly, to systems satisfying the thermodynamic limit (corresponding 
to N,V,E — > oo such that number density n — N/V and energy density 
e = E/N remain constant). However, the rapid experimental and computa- 
tional progress during the last two decades led to an increasing interest in 
extending thermodynamic concepts to 'small' systems, containing - by def- 
inition - only a very limited number of DOF [21-25]. Experiments on fi- 
nite systems include, e.g., investigations of two-dimensional Coulomb clus- 
ters in dusty plasmas [26,27], Bose-Einstein condensation in magneto-optical 
traps [28,29], and transitions in sodium chistcrs [30]. These experimental in- 
vestigations were accompanied by extensive theoretical and numerical studies 
(see, e.g., Refs. [31-33]). 

Simultaneously, interest began to focus on the question how to identify and 
classify the finite-size analogues of macroscopic PTs. Import results in this 
regard were obtained by Wales et al. [23,24], who considered necessary and 
sufficient criteria for phase coexistence in finite systems. A general classifica- 
tion scheme for smooth canonical PTs in small systems was proposed by Bor- 
rmann et al. [34]. Pursuing an approach similar to that of Yang and Lee [8,9], 
Fisher [10] and Grossmann et al. [11-13], these authors suggest to use the DOZ 
in order to characterize transitions in the canonical ensemble of small systems 
(cf. Sec. 3.2 below). Miilken et al. [35] and Alves et al. [36] compare the DOZ 
classification scheme with alternative proposals made by Gross [25,37] and by 
Janke and Kenna [38], respectively. 

Such progress notwithstanding, there still exist some open questions regard- 
ing which types of non-analytic PTs can occur in small systems. ¥or(grand- 
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)canomcal ensembles, it is well established that truly singular PTs can be 
observed in the thermodynamic limit only, corresponding to a system with for- 
mally infinite particle number N ^ oo [8-13,34,39]. In contrast to this, the 
microcanonical thermodynamic functions (TDFs) may exhibit non-analytic 
behavior even at finite N. For example, recently, Pleimling and Behringer [40] 
have found singularities in microcanonical quantities of finite three-dimensional 
spin models, which announce a continuous macroscopic PT of the infinite sys- 
tems. Additionally, as we intend to demonstrate here by means of very sim- 
ple examples, the microcanonical TDFs of a small system can also exhibit 
non-analytic microscopic PTs ^ , characterized by well-defined critical energy 
values and typically accompanied by strong variations of thermodynamic ob- 
servables (e.g. oscillations of temperature, heat capacity and pressure). From 
the physical point of view, such microscopic PTs correspond to transitions 
between different dissociation states of the system. Hence, they are important 
indicators for essential structural changes in the small system under consid- 
eration (quite analogous to singular points indicating macroscopic PTs). For 
reasons of simplicity, the discussion in the present paper will be restricted to 
one-dimensional (ID) models, but the general mechanism responsible for the 
singular and oscillatory behavior of microcanonical thermodynamic observ- 
ables works analogously in two and three spatial dimensions. Hence, it should 
in principle be possible to observe such microscopic oscillations in suitably 
designed evaporation/dissociation experiments (which must, of course, realize 
the physical preconditions of the microcanonical ensemble). 

The paper is organized as follows: Section 2 is dedicated to microscopic PTs 
in the microcanonical ensemble (MCE) of small model systems. As examples, 
we will consider isolated ID chains with Lennard- Jones (LJ) pair interactions 
and also the Takahashi gas [44]. It will be shown that these simple systems 
exhibit singular microscopic PTs, separating different microcanonical dissoci- 
ation states. Subsequently, the Takahashi gas will be used in Sec. 3 to investi- 
gate the relation between singular microscopic PTs in the MCE and smooth 
PTs in the canonical ensemble (CE) as defined by the DOZ scheme [34]. The 
paper concludes with a summary of the main results in Sec. 4. 



2 Microscopic phase transitions in the microcanonical ensemble 

Classical microcanonical thermodynamics refers to an ensemble of thermally 
isolated systems, completely described by their Hamiltonian dynamics. Due 

^ The appearance of such microcanonical singular points even in the ID case is not 
in conflict with van Hovc's theorem for the canonical ensemble [41,42], since for 
most small systems (as well as for many 'large' systems) the microcanonical and 
canonical ensembles are generally not equivalent (see, e.g., Costeniuc et al. [43]). 
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to the fact that the systems are decoupled from the environment, the energy 
E is a conserved quantity, i.e., there are no energy fluctuations in the MCE. 
In particular, as will be demonstrated in Sees. 2.3 and 2.4, the micro canonical 
TDFs may exhibit singular points even in the case of small systems. 



2.1 The microcanonical ensemble 



For the sake of simplicity we will confine ourselves to examples, where the 
thermodynamic state is completely characterized by two control parameters, 
namely, energy E and volume number V (generalizations to problems with 
additional macroscopic variables, e.g. the strength of external magnetic fields, 
are straightforward). More precisely, we will consider N identical point- like 
particles of mass m, moving in D equivalent spatial dimensions; i.e., the num- 
ber of DOF reads d — DN and the volume is given hyV — L^, where L is the 
length of the confining cube (volume interval). The deterministic dynamics of 
the system is assumed to be governed by a Hamiltonian of the standard form 



H{q, p; V) = Kip) + U{q; V) = E, (la) 

where q = {qi, ■ ■ ■ ,qd) and p = {pi, ■ ■ ■ ,Pd) are generalized coordinates and 
momenta, respectively. The kinetic energy K and the potential energy U are 
given by 



^(p) = E ' ^(^; ^) = ^p-(^) + ^box(?; V), (lb) 

where [/pair(9) represents the pair interactions of the particles, and the box 
potential is defined by 




U^oM'^V) = { . (Ic) 

otherwise. 

The primary thermodynamic potential of the MCE is the entropy 5* = S{E^ V), 
which is related to the microcanonical 'partition' function 2.m(-E, V) by [45] 



S^khlZyi, (2) 

with k denoting the Boltzmann constant. The equations of state (EOS) for 
the temperature T and pressure P are obtained by [25,45-49] 
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That is, the temperature T of the MCE is a derived quantity, which is in 
contrast to the CE (see Sec. 3.1 below), where the temperature is one of the 
adjustable external control parameters. 

Dming the past century, various different definitions for the microcanonical 
partition function Zm, or the entropy S, have been proposed and investigated. 
For classical systems as described by the Hamiltonian function (la), the two 
most commonly used definitions for Zm{E, V) read [25,45-49] 



ZM = n, (4a) 
ZM = e^ — , (4b) 

where the phase volume Q is given by {h denotes the Planck constant) 



The Heaviside unit step function Q{x), appearing in Eq. (5), is defined by 

Q{x) = for a; < and ©(x) = 1 for x > 0. The additional parameter eo 
in Eq. (4b) is a small energy constant that quantifies the thickness of a thin 
energy shell around the phase space surface defined by H{q,p, V) = E and is 
formally required to make Zu dimensionless. 

It is well-known that the definitions (4a) and (4b) may yield (almost) identical 
results [39,45,47,48] in the thermodynamic limit, i.e., if is large. However, 
for small systems they lead to essentially different TDFs. To briefly illustrate 
this, let us consider an ideal gas with N non-interacting particles, moving in 
the D-dimensional volume V. In this case, the phase volume is given by [45] 



where r(a;) is Euler's Gamma-function. Definition (4a) then yields the EOS 



2 ' TV' 
whereas one obtains from definition (4b) 



(7) 



For systems with d ^ 1 the difference in the energy equations is neghgible, but 
for small systems it becomes relevant. In particular, for d = 1 Eq. (8) yields a 
negative temperature at positive energy. Obviously, a similarly unreasonable 
result is obtained for d = 2. This indicates that the partition function (4b) is 
inappropriate for systems with low-dimensional phase space. 

More generally speaking, only Eq. (4a) reproduces correctly the well-known 
laws of thermodynamics and also yields the correct equipartition theorem for 
an arbitrary number d of DOF, whereas the Eq. (4b) leads to inconsistencies 
if d is small. This important aspect was first realized by Hertz [50,51], and, 
later on, also emphasized by Becker [45], Berdichevsky and v. Alberti [52,53] 
and Adib [54]. In particular, denoting the ensemble average with respect to 
the microcanonical probability density function f{q,p) oc 6{E — H{q,p)) by 
( ■ )mce, it can be shown [45,46] that the (equipartition) identity 



Due to these reasons, all subsequent considerations will be based on Eq. (10). ^ 
For Hamiltonians as in Eq. (la) one can still perform the momentum integra- 
tion in Eq. (5) using d-dimensional spherical coordinates, yielding 



where Oa — 27r''/^/r(d/2) denotes the surface of the d-dimensional unit sphere. 
2.2 Macroscopic vs. microscopic phase transitions 

Conventionally, macroscopic PTs are singularities (non-analyticities) of TDFs 
that arise in the thermodynamic limit, corresponding to an infinitely large 
system. The first systematic classification of macroscopic PTs goes back to 
Ehrenfest [3,55]. According to the Ehrenfest scheme, a PT is indicated by 
a non-analyticity of the Gibbs free enthalpy G {T, P, . . .), assumed to be a 
function of the temperature T, pressure P and other external control param- 
eters. The order of the PT is determined by the lowest order at which any of 

^ This is e.g. in contrast to Refs. [25,35,37] where the Boltzmann definition (4b) is 
considered. 





(9) 



holds, only if one employs the Hertz entropy definition 



S = k\nn. 



(10) 



6 



the derivatives of G (T, P, . . .) becomes non-continuous. Since the Ehrenfest 

scheme has turned out to be too narrow in many cases, it is nowadays of- 
ten preferred to merely distinguish discontinuous {first-order) and continuous 
{second- order) transitions. 

Extending this concept to small systems, any non-analyticity of the ther- 
modynamic potential as function of the external control parameters may be 
called a PT. However, to avoid ambiguities, we shall speak of microscopic PTs 
when discussing singular (non-analytic) points in the micro canonical TDFs of 
small systems. Furthermore, we will adopt the following terminology to clas- 
sify microscopic PTs in the MCE: If the primary thermodynamic potential, 
the Hertz entropy S*, is discontinuous, then we will call the PT discontinuous; 
if S is non-analytic but continuous, the microscopic PT is called continuous. 

Let us next discuss how the formal (Ehrenfest-type) order of a microscopic PT 
depends on the number of DOE. Eor systems described by Hamiltonian (la), 
the phase volume (11) is related to the admissible configuration space volume 



, 3 



{E,V)^ I dqe(E-U{q;V)) (12) 



Via 

d'^l^^{E,V) _ {2^mYI\ 

QEm -^v!7^^^^'^^- ^^^^ 

If, at a given critical energy EfXY)^ ijj{E^'V) has continuous derivatives up 
to order j, but a discontinuous (j -|- l)st derivative, then f2(£^, V) and hence 
S{E^ V) have continuous derivatives up to order (j -|- (i/2), but a discontinuity 
in the (j -|- (i/2 -|- l)st derivative. Consequently, the formal order of the PT 
equals (j-|-(i/2-|-l), i.e., the order increases with an increasing number of DOE. 
Eor the examples discussed below, the differentiability class of the admissible 
configuration space volume lo{E^ V) does not change with particle number 
(j = for Lennard- Jones chains and j — —1 for the Takahashi gas), and we 
indeed observe such an increasing order with increasing particle number (cf. 
results of Sec. 2.3.2).'' 



2.3 Singular microscopic phase transitions in Lennard- Jones chains 



To demonstrate the appearance of non-analytic microscopic PTs in the MCE, 
we consider a ID Lennard- Jones (LJ) chain, moving freely in a ID box volume 
[— L/2, L/2\. In this case, the pair potential in Hamiltonian (la) reads: 

^ In an ordinary sense, Eq. (13) is defined for even integer values d> only; how- 
ever, by employing fractional derivatives [56], its range of validity can be extended 

to odd integer values d > 0. 

^ Recently, similar results have been reported for the mean-field spherical spin 
model by Kastner and Schnetz [57]. 
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1 ^ 

where a, cr > are positive parameters. To simplify subsequent formulae, we 
will measure energy and length in units of the parameters a and tq = 2^^^a, 
where tq is the position of the minimum of t/Lj(r). With respect to these units 
the LJ potential (14) is given by 

^La(r) = ^ - (15) 

For small volumes L < 1 the LJ-force is always repulsive, whereas in the more 
interesting case of sufficiently large volumes, L > 1, the LJ-force may also 
become attracting. 



a 



12 



6 



(14) 



2.3.1 Diatomic Lennard- Jones molecule 

We start by discussing the simplest non-trivial example N = 2 and D = 1, 
where we can calculate the microcanonical TDFs exactly. ^ In this case, the 
energy E can take values Eq{L) < E < oo, where the groundstate energy is 
given by 



p.n- ■ uf \ IUl,{L), L<1, 
Eo{L) = mmH{q,p)^ i 

(g,p) 1—1, L > 1. 



(16) 



A straightforward calculation of the phase volume, based on relative and 
center-of-mass coordinates, yields 



27rm / L 



1 2L 1 _ r^E^ 
^ + ^ + LrE - —— 



\^ 11^11 lOrio 5^5 2r^ 
where the boundary values are given by 



(17a) 



^m.n A , W |^_i/6^ E<E,{L), ^ ^ 

using the convenient abbreviations X = 1 + E and Y = 1 — Vl + E. 



^ The phase volume for the more complicated three-dimensional problem was re- 
cently calculated by Umirzakov [58]. 
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Fig. 1. Microcanonical phase diagrams for the ID LJ molecules. Energies below 

the minimum energy Eq{L) (solid line) arc forbidden, (a) Diatomic LJ molecule 
{N = 2): The critical curve Ec{L) (dashed line) separates a gas- like (or dissociated) 
phase from a molecule phase, (b) Triatomic LJ molecule {N = 3): The critical curve 
Ec\ (L) (dashed line) separates a gas-like (or dissociated) phase from an intermediate 
(partially bound) phase, enclosed by the critical curves Ec\{L) and Ec2{L) (dotted 
line), and a bound molecule phase (3) below Ec2{L). 



In the case L < 1 we have rmax = L for all energies E > Eo, and, hence, the 
phase volume Q(E) is a smooth function for all permitted energy values. For 
L > 1, however, the boundary value rmax changes its energy dependence at 
E = Ec{L) in a non-analytic manner, and hence the phase volume Q{E) is not 
analytic at E = Ec{L). The critical curve Ec{L) = Ui,j{L), L > 1 separates a 
gas-like phase (dissociated state) from the molecular phase (bound state) in 
the (L, £')-parameter plane. This is illustrated in Fig. 1 and will become par- 
ticularly evident from the expansions presented in the next two paragraphs. It 
is worthwhile to stress again that the critical curve Ec{L) arises naturally due 
to the sudden change in the energy dependence of the phase volume, occurring 
when the energy E passes the critical curve. The microcanonical caloric curve 
T{E) is continuous but not differentiablc along the critical transition curve 
Ec{L), which is located in the region of an S-bend or van dcr Waals-type loop, 
respectively. Formally, this corresponds to a fourth-order transition. 



Super-critical energy values (dissociated phase). Using the result (17) 
with r2 = L, corresponding to supercritical energy values E > Ec{L) - or 
region '(!)' in Fig. 1, respectively - we can derive from Eqs. (3) the micro- 
canonical EOS, yielding 



kT 



ry I 24L(X-5£;) 



33(X-2E) 
XV3 



55 



T2 _ (g-|-X)(2LXV6_i) 

^ X4/3(l+B)l/2 



110M+^-|i 



44 



L4 



2L[X(5X-22)-h55E] 



ry , ll[X(X-5)-h5B] 
^ ^ XI73 



2L[X(5X-22)+55g] 



(18a) 
(18b) 
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where Z = 55 L'^E + IIL ^ — L Taking the high-energy hmit at constant 
volume V — L one finds 



kT PL 
lim — = 1 , lim = 2 , (19) 

corresponding to the laws for the ideal ID two-particle gas. Hence, the param- 
eter region E > Ec{L) can be identified as two-particle gas state or dissociated 
phase, respectively. 



Sub-critical energy values (bound phase). Because of Eq. (16), the 
opposite case E < Ec{L) - corresponding to region '(3)' in Fig. 1 (a) - can 
only be realized, if L > 1 holds; i.e., if the box volume is larger than the 
distance corresponding to the potential minimum tq. Again using Eqs. (17), 
but this time with r2 = Y~^/^^ we obtain 



3[8Lg(x,y)-ii/(x,y)](-i^)V3 

55 [XV3 _ yi/3 + 2LXV6(yi/3 _ (-E)i/6)] ' ^^''^> 

kT 8L^(x,y)-ii/(x,y)' ^^^^ 

where we have made use of the abbreviations 

^ = (1 + 5X)Y^'^ - (1 + 5Y)x^i^, / = (1 + 2X)y2/^ - (1 + 2y)x^/^ 

Expanding EOS (20b) near the groundstate energy Eq = —1 yields 



kT ^ -{E + 1) + Q\{E + if 



(21a) 

(21b) 



Equation (21a) indicates that, at low energy, each momentum variable as well 
as an approximately harmonic excitation of the relative coordinate carries on 
average the energy amount kT/2. Obviously, this is in agreement with the 
equipartition theorem for harmonic DOF. Equation (21b) corresponds to the 
pressure law for an ideal one-particle gas in the reduced ID volume V^s = L—1, 
reflecting the fact that, at sufficiently low energy values, the two particles form 
a bound molecule with distance r ~ 1 between each other. 



In Fig. 2 (b) the caloric curves and the pressure law are shown for different 
fixed values of L and E, respectively. When the volume is large enough, L > 
1, the caloric curves exhibit a characteristic convex region (S-bend) and, in 
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Fig. 2. Microcanonical TDFs for the ID diatomic LJ molecule (iV = 2). Energy 
is measured in units of the binding energy a. Volume V = L is measured in units 
of the parameter tq, corresponding to the minimum of the LJ potential. The mass 
unit is chosen such that m = 1. (a) Hertz entropy S' = S — Sq as function of 
the energy for L = 20 (solid line), L = 3 (dotted line), and L = 0.95 (dashed), 
where Sq = —k\n[h'^/{mrQa'^)]. Note the convex curvature of the solid curve at the 
transition energy, (b) Caloric curves for L = 20 (solid line), L = 3 (dotted line), and 
L = 0.95 (dashed). One can readily see the singularity (peak) in the S-bend region, 
occurring exactly when the critical line Ec{L) = Ui,j{L) in Fig. 1 is crossed. 

particular, also a non-differentiable point (see solid and dotted curves). These 
kinks occur when the energy passes through the critical value Ec{L). 



2.3.2 Lennard- Jones chains with N > 2 particles 

Analogous microscopic PTs do also occur in LJ molecules with larger particle 
numbers. For A?^ > 2 it is very difficult or, perhaps, even impossible to express 
the phase volume (11) in terms of closed functions. Usually, one can perform 
only one or two of the integrations analytically, and the remaining integrals 
have to be calculated numerically, using e.g. Monte-Carlo methods. We em- 
ployed the Divonne algorithm of the CuBA library [59] to calculate the phase 
volume for LJ molecules with 3 and 4 particles. We used at least 1 miUion 
sample points and partially increased the sample size up to 100 million points 
for testing. We also cross-checked the results with other deterministic and 
probabilistic integration algorithms of the CUBA library and found excellent 
agreement. 

Figure 3 shows the numerically calculated microcanonical caloric curve T{E) 
of the three-particle LJ chain (A^ = 3) and its first derivative for different 
values of L. For L = 1.9, the caloric curve appears smooth and almost linear. 
For L — 6 and = 40, T{E) is still continuous, but shows two S-bend regions 
around energies — 1 and Ec2 ~ 0. At these energies, the first derivative 

of T{E) exhibits a (negative) 'lambda peak', whereas the second derivative 
diverges. 

As in the two-particle case, the phase volume fl{E, L) can be written as an 
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integral with an integrand which is analytic for all values of E and L, and 
integration boundaries which are analytic except for certain values of E and 
L. A detailed analysis yields the two critical curves 

^ci(^) = t/u (§) + t/u (!) + t/u (L) if L > (^) , 

Ec2{L) = t/Lj {rc2{L)) + U^j (L - r.^iL)) + U^j (L) if L > 2 (^y j , 

where rc2{L) is given by a polynomial equation of degree eighteen, with rc2{L) ~ 
— 1 for L 3> 1. The resulting phase diagram is shown in Fig. 1 (b). One readily 
identifies three microscopic phases separated by the critical curves Eci{L) and 
Ec2{L). 

At high energies E > Eci{L), the system is in a gas-like, fully dissociated 
phase. The relative positions r2i = \q2 — Qil and = \q3~ q2\ of the particles 
are only restricted by the hard-core repulsive part of the interaction potential, 
but apart from this constraint the particles can move independently inside the 
remaining volume. In the high-energy limit, one finds E ^kT, corresponding 
to a quasi-ideal ID three-particle gas. 

1 /6 

For L > 2 (^Y^ and low energies Eq{L) < E < Ec2{L), the system is in a 
bound molecule phase. The relative positions r2i and are restricted by the 
interaction potential to be close to the equilibrium position. 

For L > 2 (^y) intermediate energy values Ec2(L) < E < Eci(L), the 

system is in a partially dissociated phase. One of the relative positions r2i 
and r32 is restricted to be close to its equilibrium value, whereas the other 
is only restricted by the hard-core repulsive part of the interaction potential 
and the box volume. Accordingly, one of the three particles may move rather 
independently inside the box, whereas the other two remain bound to each 
other. 

Although more complicated in detail, the calculations for the four-particle LJ 
chain are in principle the same as for A?^ = 3. We briefiy list the main results: 
For N — A (and only taking nearest-neighbour interaction into account), three 
critical lines Eci{L), Ec2{L) and Ecsi^L) divide the (L, £')-plane into four dif- 
ferent phases. At high energies or at low volumes L < Lq, Lq = 3, the system 
is in a gas-like phase. At low energies and large volume, the system is in a 
molecule phase, where the particles can only move as a whole molecule inside 
the volume. Between the gas phase and the molecule phase, there are now two 
intermediate, or partly dissociated, phases, where the LJ molecule is broken 
up into two or three parts consisting of one or two particles. For large vol- 
umes, the caloric curve T{E) shows three PTs with a continuous first, but 
discontinuous second derivative. This confirms that the formal order of the 
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microscopic PTs increases with particle number as discussed in Sec. 2.2. 



2.4 Takahashi gas 



As another example, let us consider the piecewise linear pair potential 



u{r) 



+00, |r| = 0, 

a\r\/ro, < \r\ < ro (22) 



where a and ro arc positive parameters. In the literature this model is known as 
the Takahashi gas [44] . The potential u from Eq. (22) is qualitatively similar 
to the LJ potential, but, unlike C/lj, it allows to calculate the exact TDFs 
for both MCE and CE in the case N = 2. The Takahashi gas is, therefore, 
particularly well suited for studying the differences between the two ensembles. 

To keep subsequent calculations as simple as possible, we will measure energy 
and length in units of a and ro from now on (ro now defines the range of the 
potential). With respect to these units, the potential in Eq. (22) simplifies to 



+00, |r| = 0, 
u{r) — < |r|, < |r| < 1, 

.1, kl>l, 



(23) 



and, in contrast to Eq. (16), the groundstate energy is now given by Eq[L) — 0. 
A straightforward calculation of the phase vohimc yields: 
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Fig. 3. LJ chain with AT = 3 particles, (a) The microcanonical caloric curve T{E), 
and (b) its first derivative for different values of the volume L = 1.9 (dashed line), 
L = 6 (dotted line), and L = 40 (solid line). T(E) shows two S-bend regions for 
L = 6 and L = 40. In these regions, the first derivative dT{E)/dE exhibits a lambda 
peak pointing downward (the second derivative d^T{E)/dE'^ has a pole). 
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(L^{3E-L), E>Ec{L), L < 1, 

TTTTl I 

l3L[L(E-l) + l]-l, E>E,(L), L > 1, (24) 
[e\3L-E), E<E,{L), 

where the critical energy curve is given by Ec{L) = u{L). Note that, compared 
with LJ potential from Sec. 2.3.1, we must now additionally distinguish the 
cases L < 1 and L > 1, but this is only because the piecewise linear potential is 
not differentiable at r = 1. However, analogous to the case of the LJ potential, 
the critical curve Ec{L) = u{L) separates a gas- like phase (dissociated state) 
from the molecule phase (bound state) in the (L, £')-parameter plane. This 
becomes evident from the microcanonical EOS: 



1 

kf 

P 

kf 



SE-L ■ 



3L[L(ii-l) + l]-l' 

2. I ^_ 
E ^ E-3L' 

3L-6E 



L^-SLE^ 

6L(E-l)+3 
3L[L(E-1)+1]-1' 

3 



E > E,{L), 
E > E,{L), 
E<E,{L). 

E > E,{L), 
E > E,{L), 
E<E,{L). 



L<1, 
L > 1, 

L < 1, 
L > 1, 



(25a) 



(25b) 



I, 3L-E' 

Taking the high-energy limits of Eqs. (25) at constant volume V — L one finds 



E PL 
lim tt;; = 1 , lim -— = 2 , (26) 

corresponding to the laws for the ideal ID two-particle gas. Hence, the pa- 
rameter region E > Ec{L) can be referred to as gas-like or dissociated phase, 
respectively. 

In the opposite case, E < Ec{L), a low-energy expansion of the EOS near the 
groundstate energy Eq = —1 yields 



Neglecting terms of higher order. Eq. (27) corresponds to the pressure law for 
an ideal one-particle gas. This reflects that at sufficiently low energy the two 
particles form a bound LJ-type molecule. Compared with the results of the 
preceding section, a slight difference is given by the fact that in the LJ case 
the particles have a non- vanishing distance in the groundstate, see Eq. (21b). 

In Fig. 4 (a) and (b) the caloric curves and the pressure law are shown for 
different fixed values of L and E, respectively. In particular, in Fig. 4 (a) one 
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Fig. 4. Microcanonical TDFs for the ID diatomic {N = 2) Takahashi gas. (a) Caloric 
curves for L = 20 (solid line), L = 5 (dotted), and L = 0.5 (dashed). One can see the 
jump singularity for L > 1, occurring exactly when the critical line Ec{L) = u{L) is 
crossed, (b) Isoenergetic pressure curves for E = 0.5 (solid line), E = 1.5 (dotted), 
and E = 5.0 (dashed), where /3 = (fcT)-^ The local minimum is observable only 
if E/a > 1. In the limit L ^ oo the particles behave as an ideal two-particle 
(one-particle) gas, if they are in the dissociated phase E > a (bound phase E < a). 



can see that the caloric curves exhibit a singularity, when the energy passes 
through the critical value Ec{L). More exactly, for L < 1 the caloric curves 
are continuous but not differentiable at = Ec{L), corresponding to a third- 
order PT. By contrast, in the complementary case L > 1, the caloric curves 
become discontinuous at Ec{L), corresponding to a second-order PT (which is 
a consequence of the additional singularity at r = 1 and the vanishing gradient 
of the piecewise linear potential at r > 1). 

In the next section, we are going to study the relationship between this singular 
microcanonical PT and the critical behavior of the corresponding systems in 



3 Smooth phase transitions in the canonical ensemble 



Although referring to different physical conditions, MCE and CE may yield 
(almost) identical TDFs for well-behaved systems with a large number of 
DOF d ^ oo. However, this 'equivalence' between the different statistical 
ensembles does usually not hold for systems with a small number of DOF 
[25,43,45,60,61]. Hence, one has to describe the thermodynamics of small 
systems by the statistical ensemble which actually corresponds to the given 
physical conditions. 



the CE. 
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3.1 The canonical ensemble 



Considering the CE is appropriate if the system under investigation is coupled 
to an infinite heat bath [45,62-64]. The surrounding heat bath (thermostat) 

keeps the temperature of the particles constant, but causes energy fluctuations 
5E > around the energy mean value E. Consequently, singularities in the 
thermodynamic functions may only exist in the thermodynamic limit with 
A^, £■ — > oo, such that e — E/N and n — N/V remain constant and Se — > 0; 
i.e., for a finite system with N < oo the heat bath smoothens the singularities 
due to non- vanishing fluctuations 6e. Nevertheless, it is possible to deflne and 
to classify 'smooth' PTs by virtue of the DOZ scheme [34] , discussed below. 

Given a Hamiltonian of the form (la), the canonical partition function is 
defined by 



^c(A V) = j^, I dg I dp eM-P^P; V)], (28) 

where P = {kT)~^. The external control variables are now {T,V) or {P,V), 
respectively. When Zc{P, V) is known, (mean) energy and pressure of the CE 
are obtained by differentiation 



f.--.nZc, P.-^ (29) 

where F = —kT In Zq is the free energy. For convenience, we are going to drop 
the over-bars and simply write E and P in the next section (over-bars will 
be reinstated occasionally, e.g., when comparing microcanonical and canonical 
quantities) . 



3.2 Classification of canonical phase transitions in the DOZ scheme 



According to Yang and Lec [8,9], the distribution of zeros (DOZ) of the grand- 
canonical partition function determines the PTs in the grandcanonical en- 
semble. Later on, a similar approach has been employed by Fisher [10] and 
Grossmann and Rosenbauer [11] to identify and classify PTs in the canonical 
ensemble. They considered the distribution of complex zeros of the canoni- 
cal partition function Zc{P), taken as complex function of the complex inverse 
temperature 

fi = P + IT (/3 > 0) . 

For finite systems, one can show that there are no zeros on the positive real 
axis, i.e., Q'(/3fc) ^ V/c. In the thermodynamic limit, however, certain points 
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Pc on the real /3-axis may become limiting values of the DOZ. By studying 
how the zeros condense near the /?-axis, one can characterize the PT. 

Although for small systems there are, strictly speaking, no non-analytic PTs in 
the CE, it may nevertheless be helpful to distinguish different thermodynamic 
phases [31,32]. The DOZ classification scheme of Bohrmann et al. [34] is based 
on the idea that the complex zeros (3k closest to the real /3-axis can be employed 
to estimate the DOZ behavior in the thermodynamic limit. The extrapolated 
hmiting values Pc are used to define the 'smooth' canonical PT of the finite 
system. To be more specific, one first numbers the complex zeros Pk, k — 
1,2,... of the canonical partition function Zc{(3) according to their distance 
to the real /3-axis, and then calculates the quantities 



7^ 



/32-Pi 



a 



/3c 




A; — 2, 3, . . . , 



log(r3) - log(T2) ' 

■Pi - 7n- 



(30a) 
(30b) 

(30c) 
(30d) 



These quantities characterize the DOZ near the real /9-axis. According to 
the scheme of Borrmann et al. [34], a first-order transition at temperature 
Tc = {kPc)~^ appears if a = and 7 = 0, whereas values 1 > a > correspond 
to second-order transitions, and a > 1 to even higher-order transitions in the 
original Ehrenfest classification. In order to actually observe a singular PT in 
the Ehrenfest sense, it is required that ti — > in the thermodynamic limit. 

A similar classification scheme for smooth canonical PTs in small systems, 
based on the average cumulative density of zeros of the canonical partition 
function, has been proposed by Janke, Kenna et al. [38,65]. Alves et al. [36] 
compared both approaches for a variety of model systems. In the following, 
however, we will refer to the DOZ-scheme by Borrmann et al. [34] in order to 
classify PTs in the CE of a small system. In particular, we aim to compare 
the DOZ-classification with the MCE results obtained for the Takahashi gas. 



3.3 Takahashi gas 



We first calculate the canonical TDFs for the Takahashi model from Sec. 2.4, 
corresponding to two point-like particles, N = 2, confined in a ID volume 
V — L and interacting via the (rescaled) piecewise linear pair potential from 
Eq. (23). For L < 1, we find exphcitly 
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rmx 



(e-/^^ + /3L-l), 



(31a) 



and in the opposite case, that is for L > 1, 



_ r 1 



e-^{l + l3-l3L)+l3L-l 



-{L-lf 



(31b) 



Accordingly, one obtains the following canonical EOS 

{ 3+,JL+c "'^(2,JL-3) r ^ -| 

/3[l+e/3-C'(/3L-l)] ' ^ ^ -L, 

6/3+(6-4/?L)(l-e^)+/32(L-l)[L-3+/3(L-l)1 ^ . 
/3{l+[l-/3(L-l)]^+2e'^(^L-l)} ^' -L- 

L < 1, 



1 e/3i'(/3L-l)+l' 
l+[l-^(L-l)]^+2e'^(^L-l) 



2[e/^-l+/3(L-l)] ^ L > 1, 



(32a) 
(32b) 



In the high-temperature limit, corresponding to /? — > at constant volume 
V — L, one finds 



hm BE = 1 , lim 3PL = 2 , 
i.e., the system behaves like an ideal ID two-particle gas in this limit. 

2.5 
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1 

0.5 



(33) 




►J 

CO. 



(b) 




5.0 




V'---" 1.5 


kT/a =0.5 



10 

L [ro] 



15 



20 



Fig. 5. Canonical TDFs for the ID diatomic Takahashi gas. Mean energy E and 
thermal energy kT are given in units of the binding energy a. The volume ^ = L is 

measured in units of the parameter rg, corresponding to the range of the potential, 
and the mass unit is chosen such that m = 1 holds, (a) Caloric curves for different 
fixed values of L. (b) Pressure law for different fixed temperature values kT. Note, 
that (3PL 2 foT L ^ oo, corresponding to the law of the ideal two-particle gas. 



Figure 5 shows several curves, corresponding to the thermodynamic laws (32). 
One may notice a strong local increase in the solid caloric curve of Fig. 5 (a). 
This behavior is associated with a smooth canonical PT in the DOZ scheme 
of Borrmann et al. [34]. To see this, we next determine the complex zeros of 
the complex partition function 2ic(/3 + iT, L) for L > 1 from Eq. (31). Figure 
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6 (a) shows the corresponding numerical resuhs found with Mathematica [66] 
for L = 20, by using a contour plot of the function \Zc{P + if, L)\. Since the 
zeros of the partition function are complex conjugate [34], only the zeros in 
the upper complex half-plane are shown. 

To also obtain an analytic estimate for the DOZ in the more interesting case 
L 1, we expand the partition function (31b) near L —>■ oo, and find 



l + _(e^-^-l) 



(34) 



Thus, neglecting terms of 0(L 2), the zeros of Zc{/3, L) are given by 



A; = ±l,±2, 



ProdLog 



-k, 



L-2 



exp 



L-2 



(35) 



where ProdLog [A;, 2;] is the kth solution for w in z = wexp{w). The result 
(35) is vahd for L > 1, and in Fig. 6 (b) we plotted /3fe for 1 < A; < 8. By 
comparing Figs. 6 (a) and (b) it becomes evident that for L > 20 the analytic 
estimate from Eq. (35) is in good agreement with the numerical results. 
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Fig. 6. Complex zeros of canonical distribution function Zc{/3 + ir) for the ID 
diatomic molecule with the piecewise linear pair interaction potential u{r) from Eq. 
(22). Quantities (3 and r are given in units of the inverse binding energy a~^. (a) 
Numerically determined zeros in the upper half-plane obtained with Mathematica 
for L/ro = 20. (b) Analytic results based on the asymptotic series expansion of 
Zc at L/ro see Eq. (35). As evident from these two diagrams, the analytic 

estimate agrees well with the numerically determined results. 



Using Eq. (35) we can calculate the characteristic quantities a, /5c, 7 from 
Eqs. (30), required for the DOZ classification. Table 1 shows a summary of the 
results for different values of the volume parameter L. The critical temperature 
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is given by kTc = l/(3c and Ec is obtained by inserting f3c into the energy 
equation (32a). According to the DOZ scheme, for 1 ^ L < oo we find < 
a < 1, corresponding to a canonical second-order PT. However, for L — > oo 
we observe that a — > 0, 7 — > indicating that in this hmit the transition 
converges to first order. Moreover, as evident from the last two columns in 
Tab. 1, for increasing volume L the critical temperature vanishes, Tc — > 0, 
while the corresponding energy values E{Tc) approach the value 0.5. However, 
by first taking the limit L — > 00 in Eq. (32a), we find that the corresponding 
asymptotic caloric curve is given by 

lim E = l + i3-^ = l + kT, T> 0. (36) 

L—*oo 

Inserting Tc into the rhs. of this equation and letting Tc — > 0, we obtain 
Ec — 1, which is in agreement with microcanonical result Ec{L = cxo) = 1, 
and corresponds to the dissociation energy. In particular, the latter result 
means that, in the case of a very large volume, very small energy fiuctuations 
(requiring T > 0) suffice to permanently break up the molecule. On the other 
hand, the system becomes deterministic at T = 0, and, correspondingly, the 
mean energy of the CE is then always given by the groundstatc value Eq = 0, 
representing the bound state [formally this corresponds first taking the limit 
T — > in Eq. (32a)]. The apparent convergence to 0.5 in the last column of 
Table 1 just reflects the fact that for Eq. (32a) the two hmits T — > and 
L — > 00 do not commute. 



L 


[ro] 


a 


7 


/3c [a-'] 


kTc [a] 


EiTc) [a] 


2- 


IQi 


0.054 


0.112 


3.5 


0.290 


0.922 


2- 


10^ 


0.052 


0.052 


9.0 


0.111 


0.688 


2- 


10^ 


0.034 


0.028 


14.1 


0.071 


0.620 


2- 


lO^o 


1.1 • 10-2 


9.1 • 10-3 


26.3 


0.0381 


0.56210 


2- 


IQioo 


2.4 • 10-^ 


1.1 • 10-^ 


235.7 


0.0042 


0.50643 


2- 


-j^QlOOO 


2.5 ■ 10-^ 


1.2 ■ 10-^ 


2310.3 


0.0004 


0.50065 



Table 1 

Canonical PT parameters according to the DOZ scheme [34] for the ID diatomic 

molecule with piecewise linear interaction potential from Eq. (22). For L — ^ 00 the 
numerical results suggest that a — > 0, 7 — > and kTc/ a — > 0. 



4 Summciry and discussion 

In this paper we have studied phase transition-like phenomena in small sys- 
tems, characterized by a finite number of degrees of freedom (DOF). The main 
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objective was to clarify similarities and differences that arise when consider- 
ing either the microcanonical ensemble (MCE), corresponding to a thermally 
isolated system, or the canonical ensemble (CE), corresponding to a system 
in contact with an infinite heat bath. 



4.1 Microscopic phase transitions in the microcanonical ensemble 

In Sec. 2 it was shown that the microcanonical thermodynamic functions 
(TDFs) of a small system can exhibit complex oscillatory behavior and singu- 
larities (non-analyticities). Analogous to macroscopic phase transitions (PTs), 
such non-analytic points can be interpreted as microscopic PTs in the MCE. ^ 
To illustrate the physical meaning of such microscopic PTs, we calculated the 
microcanonical TDFs for two shghtly different ID toy models (Lennard-Jones 
and Takahashi gas). For both systems one can identify critical energy curves 
Ec{L) along which the primary microcanonical thermodynamic potential - 
the Hertz entropy - is non-analytic, reflected by kinks or discontinuities in the 
TDFs. In these models the microscopic PTs separate energetically different 
dissociation states. Their number and formal order increases with increasing 
particle number N [57]. In general, our results indicate that, typically, micro- 
scopic PTs in the MCE are accom,panied by strong qualitative changes of the 
thermodynamic observables, as e.g. rapid drop-offs or oscillations of temper- 
ature and heat capacities, see Fig. 3 (a). These effects should be observable 
in suitably designed evaporation experiments, realizing the conditions of the 
MCE (similar to those of Schmidt et al. [30], but without heat bath). 

For the model systems considered in this paper, a non-analyticity in the TDFs 
is characterized by a critical energy curve Ec{L). Mathematically, such critical 
curves Ec{L) arise due to the integration over the ©-function in the definition 
of the microcanonical phase volume Q. The ©-function in Eq. (11) effectively 
constrains the range of the integration to the subset 

A^{{q^,...,qd)eR''\E- U{q^, . . . , g,; F) > O} (37) 

in the configuration space; i.e., A is the energetically permitted configura- 
tion space region. Hence, a singularity in the microcanonical TDFs may arise 
whenever A changes its geometry or 'shape' in an irregular manner during a 

^ Formally, such singularities can be softened by using an 'artificially' smoothened 
box potential. However, this would not affect experimentally observable effects, as, 
e.g.; a drop-off of the temperature (mean kinetic energy) when the energy level for 
the next dissociation step is crossed. Loosely speaking, employing a smoothened box 
potential in the thermodynamic analysis of a small system is similar to avoiding the 
thermodynamic limit when being interested in PTs of large systems. 
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small variation of the control parameters E and V. One possible origin for this 
may be a change in the topology of A (as discussed by Pettini et al. [18-20]). 
For the models considered here, however, the topology of A remains unaf- 
fected and another general mechanism is at work. To illustrate this in more 
detail, we recall the example from Sec. 2.3.1, corresponding to the two-particle 
Lennard- Jones gas. In this case, the set A can be expressed as 

^ = {{qi, (I2) e I (gi, 52) e [-L/2, L/2f A E- Ui^j{qi - 52) > o} , 

where [/lj is given by Eq. (14). The first constraint for (gi, ^2) reflects the box 
potential, whereas the second constraint arises from the interaction potential. 
For E < Ec{L), A consists of the two diagonal 'strips' 



^± = {(gi,?2) e [-L/2,L/2f A r^{E) <±(q,-q,) < r+(E)} , (38) 

where r_{E) > and r^{E) > denote the classical turning points of the 
LJ potential. The 'strips' A± are bounded by the box potential on two sides, 
whereas the other two boundaries are determined by the LJ interaction po- 
tential. For E > Ec{L), however, the regions A± become triangles, bounded 
by the box potential on two sides and by the interaction potential on the 
remaining side: 

A±^{{q,,q2)e[-L/2,L/2f A r_{E) < ± {q^ - q^)} . 

Evidently, A changes its geometry dramatically, when the energy passes through 
the critical energy Ec{L), thereby giving rise to the microscopic PT. For a suf- 
ficiently large volume L ^ ro, where ro is the range of the pair interaction, 
this non-analytic transformation of A at the dissociation energy is accompa- 
nied by a change of A^s effective dimensionality. Here, we mean by 'effective 
dimensionality' the number of orthogonal configuration space directions in 
which the set A extends comparably to the system size L. For energy values 
E <S Ec{L), the 'strips' A± defined by Eq. (38) are very narrow compared to 
the system size L, and, hence, A = A+ U A- is effectively one-dimensional. 
When the energy approaches Ec{L) from below, the (average) width of the 
two subsets A± grow very rapidly, and for E = Ec{L), they become triangles 
whose size is of the order of L^; i.e., the set A is effectively two-dimensional for 
E > Ec{L). In particular, it is this very rapid growth of ^1 - or Q, respectively 
- which leads to a negative slope of the caloric curve T{E, V) in the vicinity 
of the dissociation energy Ec{L). ^ 



This also explains why microscopic PTs appear more pronounced for larger L. 
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4-2 Comparison with the canonical ensemble 



In contrast to the microcanonical TDFs, singular PTs cannot occur in the CE 
of a finite system [8,9,11-13,41,42]. This is a consequence of the two completely 
difi^erent physical conditions, underlying MCE and CE, respectively [25,45]. In 
spite of lacking sharp transitions, it is useful to distinguish different thermody- 
namic 'phases', when considering CEs of finite systems [31,32]. Following the 
proposal of Borrmann et al. [34], we determined the distribution of complex 
zeros (DOZ) for the canonical partition function of the two-particle Takahashi 
gas (Sec. 3) and found a 'smooth' canonical second-order PT (according to the 
DOZ classification scheme), provided the volume L is large compared to the 
range of the potential, but still finite; for L ^ oo the canonical PT changes 
to first order in the DOZ scheme. 

In principle, however, there exist important differences between the smooth 
PTs in the CE and the singular microscopic PTs in the MCE: Smooth canoni- 
cal PTs can be viewed as direct finite size counterparts of macroscopic PTs. In 
particular, for particle systems as discussed in this paper the DOZ-scheme 
usually yields only one transition point. Furthermore, the DOZ-order of the 
transition point is typically close to one or two. By contrast, when considering 
the corresponding isolated system (i.e. the MCE), both the number and the 
formal order of the singular microscopic PTs increase approximately propor- 
tional to the particle number (since microscopic PTs signal each dissociation 
step separately). Such differences notwithstanding, canonical and microcanon- 
ical partition functions are linked by a Laplace transformation, which suggests 
that (not only for small systems) there might exist a direct connection between 
the canonical DOZ and the appearance of microscopic non-analyticities in mi- 
crocanonical entropy (see also [67] for a discussion of this hypothesis). 

To conclude with, if one wishes to describe small systems by means of a ther- 
modynamic approach, then the non- equivalence of the different statistical en- 
sembles [43] renders necessary to specify in advance, whether or not a heat 
bath (thermostat) is coupled to the small system under consideration. The 
discussion in the present paper has focusscd on two extreme limit cases, cor- 
responding to a vanishing heat bath (MCE) and an infinite heat bath (CE). 
In the future, it would also be interesting to study phase transition-like phe- 
nomena in small systems coupled to a finite heat bath [68-70], e.g., on the 
basis of Tsallis' generahzed statistics [71-73] . 

The authors thank W. Ebeling, D. H. E. Gross, P. Hanggi, P. Talkner and an 
anonymous referee for numerous helpful remarks and suggestions. 
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